Influence of photon-assisted tunneling on heat flow in a normal metal- 
superconductor tunnel junction 
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We have investigated theoretically the influence of an AC drive on heat transport in a hybrid 
normal metal - superconductor tunnel junction in the photon-assisted tunneling regime. We find 
that the useful heat flux out from the normal metal is always reduced as compared to its magnitude 
under the static and quasi-static drive conditions. Our results are useful to predict the operative 
conditions of AC driven superconducting electron refrigerators. 
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I. INTRODUCTION 

Photon-assisted tunneling (PAT) has been dis- 
cussed in literature for almost half a century by 
2 . 3.4. 5. 6. 7. 8. 9 . 10. 11. 12. 13. 14. 15. 16. 17. 18 This phe- 

nomenon arises when a relatively high frequency field 
is applied across a tunnel junction whose DC current- 
voltage characteristics are highly non-linear. The 
radiation field is, however, slow enough to guarantee 
adiabatic evolution of the energy levels of the electrons. 
A typical system to observe PAT is a SIS tunnel junction, 
with superconducting (S) leads and a tunnel barrier (I) in 
between. Even though the system as such is a Josephson 
junction for Cooper pairs, PAT deals with the influence 
of the radiation on quasiparticle tunneling. Our system 
of interest here is a NIS tunnel junction, where one of 
the conductors is a normal metal (N). Such junctions 
exhibit highly nonlinear current-voltage characteristics 
at low temperatures, and normally the current is due 
to quasiparticles only. NIS-junctions are known to have 
peculiar heat transport properties under the applica- 
tion of a DC bias voltagei^'^i^'SS^'^s^^^^, or 
(" quasi-static" ) AC radiation of relatively low frequency 
in form of either periodic or stochastic driv o'^°i'^^'^^ . 
Specifically, it is possible to find operation regimes where 
the normal metal is refrigerated and the superconductor 
is overheated, and in some special situations the opposite 
can occur as well. The question remains whether and 
under what conditions the relatively high frequency 
radiation responsible for PAT would either enhance 
or suppress the thermal transport in the NIS-system. 
In this paper we show that the influence of PAT, as 
compared to static and quasi-static AC drive conditions, 
is to decrease the refrigeration of the normal conductor, 
and also to change, usually to increase, the magnitude of 
heat dissipation in the superconductor. Although these 
results are somewhat unfortunate for high-frequency 
applications of NIS-junctions, they are, however, useful 
in finding operating conditions, for instance, for AC 
driven electronic refrigerators22i2i. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the theoretical framework together with the dis- 




FIG. 1: The system under investigation is composed of a su- 
perconductor (S) tunnel-coupled to a normal metal (N) layer 
through an insulating barrier (I) of resistance Rt. The super- 
conductor is AC voltage biased with X4c cos(cjoi), while the N 
electrode is biased with a static voltage U. Both electrodes 
are assumed to be in thermal equilibrium. 



cussion of the conditions of its validity. In particular, in 
Sec. Ill Al we present our analytical results for the heat 
and charge currents. In Sec. Illll we show and discuss the 
results. Finally, our conclusion are drawn in Sec. IIVI 



II. MODEL AND FORMALISM 

The system under investigation consists of supercon- 
ducting (S) and a normal (N) electrode tunnel- coupled 
through an insulating barrier (I) of large resistance Rt. 
An AC voltage bias ifs, of frequency vq = wq/Stt and am- 
plitude Vac, is applied to the S electrode, while a static 
voltage (fN = U is applied to the N contact. The total 
voltage across the junction is (pN — ips = U — Vac cos ujot. 
One could, of course, consider both the AC and DC volt- 
age to be applied to the normal lead, instead. However, 
we choose the setup as shown in Fig.[l]to directly demon- 
strate equivalence of the two connections when one of the 
leads is in the superconducting state. 

In the tunnelling limit with large resistance Rt the cur- 
rents through the contact are small. If the AC frequency 
is small compared to the superconducting gap, ujq <^ A, 
the deviation from equilibrium in each lead is negligible. 
In particular, the equilibrium is preserved with respect 
to the superconducting chemical potential fis in the S 
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electrode (which has dimensions much bigger than the 
branch-imbalance relaxation length). This leads to the 
standard assumptioi>i^^'i^ 



hdx 
2 dt 



(1) 



where x is the order parameter phase. 

In the case of equilibrium described by Eq. ([T|), the 
order parameter has the form 



A(r, t) = Ao(r) exp /i, dt'^ 



(2) 



It is convenient to start with the Bogoliubov-de Gennes 
equation (BdGE) for the eigen-functions of the system, 



ih 
ih 



du 
dt 
dv 
dt 



(Ho + eips)u + Av , 
-(i?o + eips)v + A*u 



where Hq is the normal-state Hamiltonian. Solutions to 
the BdGE have the form 

u(r,i) = yo(r)e-'^*/'^+('/'^)/'^^''*' , (3) 
vir,t) = „o(r)e~'-^*/'^-(*/'*)-/'*^='''*' , (4) 

where uo(r) and vo{r) satisfy the BdGE in the absence 
of the applied potential {(ps — 0) 

Euo = HqUo + A{r)vo , 
Eva = -HaVf) + A*{r)uo . 

Using the standard approach we define the retarded 
(R) and advanced (A) Green function, which can be writ- 
ten as a matrix in Nambu space: 



qR(A) pR(A) 



where F-^^'^'> refer to the anomalous Gorkov function. 
Since these functions are statistical averages of the par- 
ticle field operators which can be decomposed into the 
wave functions Eqs. ^ and ([5]), the retarded and ad- 
vanced Green functions with help of Eqs. ^ and (jH) 
take the form 

where G^^-^^'°, F^(^)'0 refer to ips = 0. 

If the AC voltage is applied to the superconductor 



Ms 



-el4c cos • 



Using the identity 



n— + 00 



-inujot 



(5) 



(6) 



where a = eVac/fu^o, and J„ is the nth order Bessel 
function. The Green functions in the frequency repre- 
sentation take the form 



G 



R(A) 



-hu) 



Jn{a)Jm(.a)G\ 



R{A)fl 

e — nhuJo,e — huj — mhujQ 



If (fs = 0, there is no time-dependence and CJe^^ea — 
27rM(ei — e2)Gej . Here and in what follows one frequency 
subscript refers to the static Green function. 

The semi-classical Green functions are defined as the 
Green functions in the momentum representation inte- 
grated over the energy variable = /2m — Ep, 



-\-oo 



Gei,e2(P,P - k) 



Under the AC drive we thus have 



MA) 
J e.e — huj 



n,k 

E 



2n5{uo - fc^o) J„(«) Jfc-n(a)/i5;:'L„- (8) 



For the Keldysh functions we use the standard 



representatioi 



.33.34 



in terms of /i and /2 which are the 



components of the distribution function respectively odd 
and even in (e,p). In Nambu space. 



+ 00 



de' 

2TTh 



[gt.Ah 



+ hf2 



(/l,ei,c', +T3/2,ei,e',)5, 



where 



T-3 



1 

-1 



In what follows we omit the integration limits if the in- 
tegration is extended over the infinite range. With Eqs. 
([7]), ([5]) the Keldysh Green functions take the form 



ye.e — huj 



rK 
J e,e- 



2j27r(5(cj - kuJo)Jn{a)Jn-kia) 

71, k 

X - ] 

[je—nhuJo iJe — nhujQi 

X [/i(e - nhujo) -I- /2(e - nhujo)] , 
^2Tr6{uj - kuJo)Jn{a)Jk-n{a) 

n,m 

~{fl,e-nhu)o + f2,e-nhuio) f t-nhuo\ > 



(9) 



where the distributions /i and /2 in the superconductor 
refer to the state with ^ps — 0. Equations for and 
are obtained from the corresponding equations for 
gR{A) ^^^^ gK substituting g ^ g, loq ^ —i^o, and 

These solutions describe a quasi-equilibrium state with 
a time-dependent chemical potential Eq. ([5]) . In the limit 
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LjQ 0, using Eq. we have gf^"^'^^ gfl^f \ 

-R{A,K) _^ -R{A,K) ^Yiich agrees with the constant- 
voltage hmit^-. 

Here we need an obvious remark. It can be shown (see 
Appendix \K\i that for /i^ satisfying Eq. ([5]) 

^27rJ(tj - kuJo)Jn{a)Jn-k{oL)<^{e - nhujo) 



whereas the heat current flowing into the electrode i is 



$(e + Ais(i))e*"*di 



(10) 



for any function <l>(e) which has no singularities. For a 
function with singularities (or large higher-order deriva- 
tives) at certain e, Eq. pU]) holds only in the limit c^o 0. 
If the density of states gf'^"^^ , f^'"^^ , and the distribution 
function were smooth functions, the quasi-static limit 
would hold for any loq; in this case all the quantities 
would simply adiabatically depend on the AC potential 
Vac- However, due to a strong singularity at e = A of the 
density of states and/or a sharp dependence of the dis- 
tribution function for low temperatures, the quasi-static 
picture breaks down for a finite ojq, determined by the 
smallest scale of the non-linearity. As a result both the 
tunnel current and the heat flux for a flnite frequency de- 
viate strongly from the quasi-static behavior. In practice, 
the reservoirs are not perfect. In particular relaxation in 
the superconductor is still an open issue. This is a ques- 
tion that needs to be addressed separately. The present 
treatment gives the answers in the case of ideal reservoirs. 

Consider the self-consistency equation for the order pa- 
rameter. Since /2 = for (ps = 0, the self-consistency 
equation A{uj) = (A/4) / //^_^(^ de takes the form 

A(cj) = Jfc(2a)27r(5(cj - /ccjo) (H) 

k 

which is the Fourier transform of Eq. Q where 

Ao(r) = (A/4) I de [/f - /^^] /i(,) 

with /i(e) = tanh(e/2T) is the order parameter for zero 
AC field. This imphes that the Eqs. are consis- 

tent. In obtaining Eq. pT|) we use 

n— + 00 

J2 Jk+n{t)Jn{z) = Jk{t - Z) . (12) 



A. Charge and energy currents 

For two tunnel-coupled electrodes, the charge current 
that flows into the electrode i is given by^ 



m 



-Ti- 



(13) 



Q(^) 



-Tr 



d_ 

di 



d_ 

dt' 



2ieipi 



-T3 



(14) 

Here i^i is the normal-state density of states in the elec- 
trode i, and ipi are its volume and electric potential. 
The collision integral in the electrode that appears 

in Eqs. (fT3)) . ([M]) contains contribution due to tunnelling 
from neighboring electrode and the electron-phonon con- 
tribution, = + The electron-electron in- 
teractions drop out from the energy current because of 
the energy conservation. The energy flow into the elec- 
trode can thus be separated into two parts. One part 
containing I^_pi^ is the energy exchange with the heat 
bath (phonons). The other part contains the tunnel con- 
tribution //^ and is the energy current into the electrode 
through the tunnel contact. The tunnel collision integral 
for the electrode 1 in contact with an electrode 2 has the 



if (1) = [5^(2) o .9^(1) - 5^(1) o g^{2) 

+5^(2) 05^(1) -g^(l)o. 9^(2)] . (15) 

Here the arguments i = 1 or 2 refer to the electrodes S 
or N. The symbol o is the convolution over the internal 
variables 

A(l)oB(2)= J A{l;h,t')B{2-t',t2)dt' . 
The factor 



parameterizes the tunnelling strength between the elec- 
trodes, Rt being the tunnel resistance. Since in the nor- 
mal state 

9§{e;h,h) = 2[/f(e)f3-H/2^(e)]<5(ti-t2) , 
the collision integral in the superconductor is 

+2[f3/f + mgi - 25f [/f r3 + m . (16) 

The even and odd components of the distribution func- 
tion correspond to the absence of the AC potential. They 
are, respectively. 



/f(e) = -n, + (l-n_,) 



= ?iAr(e + eU) — nN{e — eU) , 



(17) 



= l-nN{e + eU)-nN{e-eU) , (18) 
for the normal lead, and 

/f(e) = -ne + (l-n-e) = 0, (19) 

/f (e) = -n, n_, = 1 - 2ns (e) = tanh , (20) 
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for the superconducting lead. Here nM{() and ns(e) are 
the Fermi functions with temperatures T/v and Tg, re- 
spectively. The distributions in the superconductor thus 
correspond to the zero-potential state. 



As far as the NIS junction is concerned, consider first 
the charge current into the superconductor defined by Eq. 
(fTS)) . The tunnel current in the frequency representation 
becomes 



Isi.^) = -^Jdeg^Y, Jn[a)Jn-k{a) {2n5{Lo - fcu;o)[/f (e) + fl{e) - /f (e + nhioo) - f^ie + nhwo)] 

n,k 
n,k 

X {2TrS{uj — kujQ)[nN{e — eU -f nhw^) — ns{e)] -)- 2'k6{uj -f kiUQ)[nN{^ — eU — nhwo) — ns(e)]} . (21) 

In Eq. (|2T|) we use the relation g^'"'^^ = —gf'"^^ for static functions, and denote g^ = {gf — gt)l^ the ratio of the 
superconducting density of states to that in the normal state, g^ — Ns{e)/N]\[. The loq component of the current is 

cos(a;ot) 



eRt 
cos(cjoO 



[deg^Y, Jn{a)Jn^i{a) [/f (e) - f^{e + nfi^o)] 

•' n 

I / dege 2, Jnict)Jn-i{ct) [?T-Jv(e + ell + nftwo) + ^^{e — eU + nftwo) ^ 2rt5(e)] 

■t J „ 



(22) 



The time averaged current takes the form 

~ls = dtg,Y.'^l{a){]^{e)-j^{e + nhw^)\ 



X [rijv(e — ell + nhujo) — ?^Ar(e + eU + nhuja)]. (23) 



The terms with /2 drop out of Eq. ([22|l due to the prop- 
erty of the Bessel functions 



J-r 



i-lTJn 



(24) 



Note that if we set uq — the average current assumes 
the zero- ac- voltage form 



(0) 



1 



2eRt 



dege[nN{e — eU) — nj^{e + eU)] 



Indeed, when taking the limit — > one should keep in 
mind that the sum ^n(ct) = 1 [which is a consequence 
of a more general relation Eq. converges at rt ~ 

a — eVac/huJo- Therefore, hnujQ ~ eVac in Eq. (|2ip : thus 
one has to put eVac — * to neglect nhujQ. However, the 
true static expression is defined for luq — but Vac 7^ 0. 
According to Eq. pi7| . it has the energy-shifted density 
of states ge±eVaa ^^'^ the distribution functions /f (e ± 
eVac), fii^ i eVac)- This static limit (i.e., t^o = and 
Kc 7^ 0) is indeed obtained from Eq. (pij) using Eq. pH)) . 
Making shifts of the integration variable we find 



2eRt 



(25) 



where fi = and 

f2 = ?^Ar(e + eC/ - eKc) -nN{e- eU + eVac) 

which corresponds to the total voltage U — Vac, according 
to Eq. (HZl). 

The heat current that flows into the superconducting 
lead can be calculated with help of Eqs. ^ and We 
find in the frequency representation 



d d 2ieips 
dt^dt'^ h 



it,t') 



(26) 



-huj 



2ieips 



d__d_ 
dt~df 



(27) 



E—hu> 



= 47rz^(5(u; - kLUo)Jn{a)Jn-k{c^){e - nhwo) 



7e—n 



and similarly for g with the substitutions g g, ips 
~ips, — *■ —Wo, and /2 — > — /2- Here the distribution 
functions again correspond to zero AC potential. 

Shifting the energy variable under the integral, the av- 
erage heat current into the superconductor becomes 

X [f^ (e) + nniuo)]de. (28) 
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The heat current Eq. ((28)) is even in ujq. For ^0 = 
Eq. ([^5]) formally goes over into 

QT = ^J^3. [/f(e)-/f(e)] de 

with /f and /f from Eqs. This is the zero- 

ac-voUage result. 

The static expression is obtained from Eqs. pO]). ([26]), 
and dSll) 

Qr*" = ^ / [/f (e) - /f de (29) 

where f^stauc^^-^ corresponds to the total voltage U—Vac, 
^iv static = l-nNie + eU-eVac)-nN{e-eU + eVac) ■ 

This should be compared to Eq. (fT5)) . 

It is also interesting to define the quasi-static regime, 
which is obtained by averaging the static heat flux 
over the sinusoidal voltage cycle with V^c — > Vac cos wo^- 
It does not coincide with the static expression due to the 
voltage oscillations. This quasi-static regime corresponds 
to the classical limit occurring at small frequencies, for 
which the photon energy fkuo is much smaller than the 
energy scale over which the non-linearity of the I-V curve 
occursi^ii^. In the system under investigation such en- 
ergy scale is set by the temperature or by the width of the 
superconducting DOS peak near the gap energy, which 
smear the sudden current onset occurring at the super- 
conductor gap. As it will be confirmed in Sec. IIIIl the 
quasi-static regime occurs for hujo <C /cb?". 

We now consider the heat current flowing out of the 
normal electrode: 

QT = Qs-{^n-^s)Is, (30) 

where Is is the tunnel charge current reported in Eq. 
(PT|) . Note that the heat extracted from the normal elec- 
trode and the heat entering the superconducting lead dif- 
fer by the energy absorbed at the NIS interface where the 
potential drops by — fs- The time- average heat cur- 
rent is 

Q°i^' ^Qs- Uls - Pac (31) 



where Pac = Kic cos{u!ot)Is ujo{t) is the average AC power 
absorbed at the NIS contact, 

Pac = / 3ey'>/n(a)J«-i(a) 

2eRt J ^ 

X [/f(e)-/f'(e + "^o)] rfe. (32) 

Note that P^c is finite both in the static case (wq = 0) 
and in the quasi-static regime (small but finite loq). 



III. RESULTS AND DISCUSSION 

We shall now discuss how the heat current depends on 
the various parameters of the system. This can be done 
by numerically evaluating the expressions given in the 
previous section. In the following we shall assume param- 
eters typical of aluminum (Al) as S material, with critical 
temperature Tc=1.19 K. We assume the superconduct- 
ing gap to follow the BCS relation Aq = UGAksTc and 
choose 



Nsie) = |Re[(e + zr)/V(e + »r)2_A2]| , 

where F is a smearing parameter which accounts for 
quasiparticle states within the gap^^^i^. We will use 
F = 10"''' Ao, as experimentally verified in Ref. [s^- Fi- 
nally, we shall always assume the N and S electrodes to 
be at the same temperature, i.e., T5 = = T. 

For the sake of definiteness, let us first consider the 
situation in which no bias is applied to the normal is- 
land {U — 0). In Figs, mja) and (b) the time-averaged 
heat current entering the S electrode (Qs), Eq. ([28]). is 
plotted as a function of the AC voltage at, respectively, 
large (T = O-SAo/fce) and small (T = 0.03Ao//cb) tem- 
peratures. The various curves refer to different values 
of frequency = luo/2tt, and calculations were per- 
formed up to vq — 40 GHz, corresponding roughly to the 
value of the superconducting gap {Aq/H ~ 43.7 GHz). 
We note that a driving frequency corresponding to 2Ao 
would lead to breaking up of the Cooper pairs. For a 
comparison we have included static (i.e., vq — 0), Eq. 
((29)) . and quasi-static regimes. For large temperatures 
[i.e., T = 0.3Ao/A:b, Fig. OJa)], the heat current Qg is a 
monotonic, nearly parabolic, function of 14c for all values 
of frequency. The first observation is that the static heat 
current is always larger than the heat current at finite 
frequency. On the one hand, it is obvious that the quasi- 
static curve is below the static one, the former being just 
an average over a cycle of the static limit (see Sec. Ill A| ). 

On the other hand, the photon-assisted heat current is 
always larger than quasi-static characteristic. To be more 
precise, the heat current monotonically decreases by de- 
creasing frequency, eventually reaching the quasi-static 
limit for small enough vq (note that the curves relative 
to z/q = 1, 5, 10 GHz arc indistinguishable from the quasi- 
static one). This means that photon-assisted processes 
give rise to an enhancement of the heat current enter- 
ing S with respect to the quasi-static situation, though 
remaining well below static values. Such enhancement re- 
flects the increase in current due to photon-assisted pro- 
cessesi^: electrons are excited to higher energy states, 
thus favoring tunneling above the gap. Of course, such 
mechanism is more effective for small temperatures. In 
such a case [i.e., T — 0.03Ao/fcB, Fig. [2Ib)], indeed, 
static and quasi-static curves present an activation-like 
behavior, with a switching voltage of T4c — 0.9Ao/e 
and Vac — 1-OAo/e, respectively, and thereby increas- 
ing almost linearly. Photon-assisted Qg increases more 
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FIG. 2: (color online) Normalized time-average heat current 
into the S electrode as a function of the amplitude of 
the AC voltage Vac for different values of z/q at (a) high 
(T = 0.3Ao/ifcB) and (b) low (T = 0.03Ao/fcB) temperature. 
The static case and the quasi-static limit are plotted for com- 
parison. Note that in (a) the curves relative to =10, 5 and 
1 GIfz coincide with the quasi-static one. 
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FIG. 3: (color online) Normalized time-averaged heat current 
out of the N metal Q^"* ^ ^ function of the amplitude of the 
AC voltage 14c for different values of frequency uo at (a) high 
(T = 0.3Ao/fcB) and (b) low (T = 0.03Ao/fcB) temperature. 
The static case and the quasi-static limit are plotted for com- 
parison. At high temperature the quasi-static limit is a good 
approximation at frequencies as low as 1 GHz. 



smoothly as compared with the quasi-static case, which 
is approached by decreasing v^. 

We now consider the heat current extracted from the 
N electrode , which differs from Qg by the AC 
power Pac (for C/ = 0) absorbed by the NIS contact [see 
Eqs. dSIl), (ESI)]- In Figs. [3i;a) and (b) we plot Q^* as 
a function of T4c for several frequencies for large and 
small temperature, respectively. The effect of Pac on 
the behavior of the heat current is very strong, giving 
rise to a maximum located around T4c = Ao/e, and to 
a sign change. By increasing Vac the heat flow out of 
N increases up to the maximum and thereafter rapidly 
decreases to negative values (heat current enters the N 
electrode). For reasons given above, the maximum quasi- 
static heat current is always smaller than the maximum 
of the static one. Another effect of Pac is that, in this 



case, the photon-assisted heat current is smaller than the 
quasi-static characteristic. In particular, the heat current 
monotonically increases by decreasing frequency, eventu- 
ally reaching the quasi-static limit for small enough vq. 
Moreover, by increasing the frequency the maximum of 
Qn moves toward smaller values of Kc • While at large 
temperatures Qn remains positive (implying heat ex- 
traction from the N electrode) also for frequencies slightly 
above 40 GHz [see Fig. [3] (a)], at low temperatures the 

minimum frequency for positive is drastically re- 

duced (by about one order of magnitude) [see Fig. [3] 
(b)]. This clearly proves that photon-assisted tunneling 
is detrimental as far as heat extraction from the N elec- 
trode is concerned. Analogously to what happens for the 
charge current^, the approach to the quasi-static limit 
depends on temperature. Indeed, as already mentioned 
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FIG. 4: (color online) The same as for Fig. [3](b) for frequen- 
cies in a smaller range. At low temperature the quasi-static 
limit is a good approximation for frequencies as low as 0.1 
GHz. 



in Sec. Ill Al the quasi-static regime occurs at hio^ <^ k-gT. 
The curve relative to 1 GHz differs, with respect to the 
quasi-static one at its maximum, by less than 0.1% at 
T = 0.3Ao/A:b, and by about 50% at T = 0.03Ao//cb, 
where k^T ~ hujQ. Figure [3] shows the time- average 
Qn versus Vac at low temperature calculated for fre- 
quencies in a smaller range. As it can be clearly seen, 
the quasi-static curve appears to be a good approxima- 
tion for i^o = 0.1 GHz. 

It is now interesting to analyze the behavior of dy- 
namic heat transport in the NIS junction for fixed am- 
plitude of the AC voltage by plotting the heat currents 
as a function of the period of oscillations tq = l/j^o- This 

is shown in Figs. [S] and [5] for Qg and Q^"*, respectively. 
Here we set C/ = 0. Both for large [T = O.SAo/fce, 
see Fig. Ela)] and small [T = 0.03Ao/fcB, see Fig. E^b)] 
temperatures the heat current Qg presents an overall de- 
crease with To, for all values of Kic- At small temper- 
atures, however, the heat current shows an additional 
structure consisting of superimposed oscillations due to 
photon-assisted processes, which tend to disappear for 
large values of tq (small frequencies), i.e., approaching 
the quasi-static limit. Notably, the relative maxima turn 
out to be equally spaced by the time scale related to the 
superconducting gap, ta = h/Ao. In addition we found 
that, at even lower temperatures, also the relative max- 
ima are equally spaced by ta- 

Though presenting an overall enhancement with tq, the 
behavior of the time-average heat current extracted from 
the N electrode is qualitatively similar to that of 

Qg (see Fig. [6l note that the vertical axis is linear in this 
case). Note that for large temperatures Qn remains 
positive for most of the frequency range considered, even 
for Vac = Ao/e [see Fig. HJa)]. For small temperatures 
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FIG. 5: (color online) Normalized time-average heat current 
into the S electrode Qg as a function of the period of the 
oscillations tq = l/uo for various values of Vac at (a) T = 
0.3Ao/fcB and (b) T — 0.03Ao/fcB- The distance between the 
relative maxima at low temperature [panel (b)] is ta = h/Ao. 



[see Fig.m^b)], however, the heat current is negative over 
nearly the whole time range. The additional structure, 
in this case, shows equal spacing (of magnitude ta) be- 
tween the relative minimums, since these correspond to 
maximum heat absorption by S [maxima in Fig. \Mjo)] ■ 

We now turn to the effect of a finite DC voltage U 
combined with an AC modulation on the heat current 
exiting the N electrode. In Figs.[7fa) and (b), the time- 
average Q™* at large temperatures (T = O.SAo/fce) is 
plotted as a function of U for several values of frequency 
at Vac — 0.6Ao/e and Vac = 0.9Ao/e, respectively. Fig- 
ure Efa) shows that is nearly constant having a 
weak maximum around U ~ 0.3 Ao/e almost indepen- 
dently of the frequency, and rapidly decreasing there- 
after. By inspecting Fig. [3ja) it clearly appears that, at 

T = 0.3Ao/fcB, Qn is maximized around T4c ~ 0.9Ao/e, 
so it seems that a finite value of U just adds to the AC 
voltage making the heat current to move along the volt- 
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FIG. 6: (color online) Normalized time-average heat current 
out of the N electrode Qn"' as a function of the period of 
the oscillations ro = l/vo for several values of Vkc at (a) T — 
0.3Ao/fcB and (b) T — 0.03Ao/fcB. At low temperature [panel 
(b)] the distance between relative minimums turn out to be 
TA = h/Ao. In (b) the curve relative to Vkc = 0.2Ao/e is 
consistent, in the scale of the figure, with zero. 



age characteristic similarly to the pure AC case. Further- 
more, we note that the addition of a static DC potential 
to an AC modulation is not able to recover the maxi- 
mum value the heat current can achieve with only the 
AC voltage biasing. A confirmation of this is given in 
the plots displayed in Fig. [7jb) which are relative to a 
value of Vac = 0.9Ao/e. For such an AC voltage biasing 

Qjyf does not present a constant part, and the addition 
of U turns out to only suppress the time-average heat 
current. Moreover, an increase of frequency i^q causes a 
reduction of Qj^ , even to negative values. 

We finally plot in Fig. [8] the maximum value of Q™*; 
obtained by spanning over V^c , as a function of T for sev- 
eral values of frequency. For every the time-average 
Qjyf is a bell-shaped function presenting a maximum 
around T ~ 0.25Ao/fcB (similarly to what happens in 
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FIG. 7; (color online) Normalized time-average heat current 
exiting the N electrode Q^"' as a function of the DC voltage 
U for various values of vo at (a) Vkc = 0.6Ao/e and (b) Vac = 
0.9Ao/e. In (a) and (b) the calculations were performed by 
setting the temperature at T = 0.3Ao/A:b. 



the static^^ as well as in the quasi-static limit), which 
is gradually suppressed upon enhancing the frequency. 
By increasing the frequency the curves slightly shrink, 
thus reducing the temperature interval of positive heat 
current. Moreover, the position of the maxima tends to 
move to higher temperatures for intermediate frequencies 
(i.e., for vq below ~ 20 GHZ), while they tend to move 
to lower temperatures in the higher range of frequencies 
(see for example the curve corresponding to i/q — 40 GHz 
in Fig. il). 



IV. CONCLUSIONS 

In this paper we have calculated the heat currents 
in a normal/superconductor tunnel junction driven by 
an oscillating bias voltage in the photon-assisted tun- 
neling regime. We have found that the maximum heat 
extracted from the normal electrode decreases with in- 
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FIG. 8: (color online) Normalized time-average heat current 
out of the N electrode Q^" , maximized with respect to Vkc, 
as a function of temperature T at i7 = for different values 
of frequency vq. The static case and the quasi-static limit are 
plotted for comparison. 

creasing driving frequency. We checked that for small 
frequencies {huj^ <^ k^T) the photon-assisted heat cur- 
rent approaches the quasi-static limit, the latter being 
obtained by averaging the static heat current over a si- 
nusoidal voltage cycle (relevant for sub-GHz frequencies). 
The suppression of the heat current by photon-assisted 
processes can be imputed to the AC power, dissipated 
at the tunnel contact, which is enhanced in the quantum 
regime with respect to the quasi-static limit. On the 
contrary, the heat current entering the superconducting 
electrode slightly increases with increasing frequency. We 
also found that, for small temperatures, the heat current 
as a function of the inverse of frequency presents an ad- 
ditional structure consisting of superimposed oscillations 
with a period corresponding to the time scale derived 
from the superconducting gap, ta = h/ Aq. 

We want finally to briefly comment onto some implica- 
tions of the above results for practically realizeable sys- 
tems. We refer, for instance, to AC-driven NIS electron 
refrigerators operating in the regime of Coulomb block- 
ade which were theoretically investigated in Rcf."^", and 
experimentally demonstrated in Refi^. More in partic- 
ular, it was shown in Ref.'^" that both the heat current 



flowing out the N island and the minimum achievable 
electron temperature depend on the frequency of the 
gate voltage as well as on the bath temperature. Our 
results may thus suggest the proper operating conditions 
in terms of frequencies and bath temperatures in order 
for photon- assisted tunneling not to suppress the heat 
current in these systems. In other words, they allow to 
predict a suitable range of parameters which keep the 
system in the quasi-static limit. 
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APPENDIX A: STATIC LIMIT 

To prove Eq. (fTU]) we use that <I>(x) is analytic, there- 
fore 

k=0 

We next perform the inverse Fourier transform of the 
right-hand side of Eq. PH)) 

$(e + fisit)) = Jn{a)Jn~k{c^)H^ - nnwo)e-''''^°' 

71, k 

n 

= ei fo <'Vaa cosc^ot' dt' ^ $ U - ih^ \ J n{a) ''''^ 

n ^ ^ 
— gTi lo COS u:ot' dt' ^ ( £_ j/j^ I ^-i Jo eVc,a COS uiot' dt' 

V dt) 

which indeed is $(e — eVac coso^o^). Here we use the ex- 
pansion Eq. ^ twice. 



* Electronic address: kopnin@boojum .hut.fi| 
t Electronic address: giazotto@sns.it 

^ A. H. Dayem and R. J. Martin, Phys. Rev. Lett. 8, 246 
(1962). 

^ C. F. Cook and G. E. Everett, Phys. Rev. 159, 374 (1967). 
^ E. Lax and F. L. Vernon, Phys. Rev. Lett. 14, 256 (1965). 

* C. A. Hamilton and S. Shapiro, Phys. Rev. B 2, 4494 
(1970). 

^ T. Kommers and J. Clarke, Phys. Rev. Lett. 38, 1091 



(1977). 

® A. A. Kozhevnikov, R. J. Schoelkopf, and D. E. Prober, 

Phys. Rev. Lett. 84, 3398 (2000). 
^ A. Vaknin and Z. Ovadyahu, Europhys. Lett. 47, 615 

(1999). 

* K. W. Yu, Phys. Rev. B 29, 181 (1984). 

^ S. P. Kashinje and P. Wyder, J. Phys. C: Solid State Phys. 

19, 3193 (1986). 
° F. Habbal, W. C. Danchi, and M. Tinkham, Appl. Phys. 



10 



Lett. 42, 296 (1983). 
" J. E. Mooij and T. M. Klapwijk, Phys. Rev. B 27, 3054 
(1983). 

^2 B. Leone, J. R. Gao, T. M. Klapwijk, B. D. Jackson, W. 
M. Laauwen, and G. de Lange, Appl. Phys. Lett. 78, 1616 
(2001). 

13 Y. Uzawa and Z. Wang, Phys. Rev. Lett. 95, 017002 

(2005). 

1* P. K. Tien, and J. P. Gordon, Phys. Rev. 129, 647 (1963). 
1^ J. R. Tucker and M. J. Feldman, Rev. Mod. Phys. 57, 1055 
(1985). 

J. R. Tucker, IEEE J. Quantum Electron. 15, 1234 (1979). 
" J. N. Sweet and G. I. Rochlin, Phys. Rev. B 2, 656 (1970). 
1* U. Zimmermann and K. Kreck, Z. Phys. B 101, 555 (1996). 
1® Sec F. Giazotto, T. T. Hoikkila, A. Luukanen, A. M. Savin, 

and J. P. Pckola, Rev. Mod. Phys. 78, 217 (2006), and 

references therein. 
^° A. Bardas and D. Averin, Phys. Rev. B 52, 12873 (1995). 
^1 M. M. Leivo, J. P. Pekola, and D. V. Averin, Appl. Phys. 

Lett. 68, 1996 (1996). 
^2 M. Nahum, T. M. Eiles, and J. M. Martinis, Appl. Phys. 

Lett. 65, 3123 (1994). 
^3 A. M. Clark, A. Wilhams, S. T. Ruggiero, M. L. van den 

Berg, and J. N. UUom, Appl. Phys. Lett. 84, 625 (2004). 
2* A. M. Clark, N. A. Miller, A. Williams, S. T. Ruggiero, G. 

C. Hilton, L. R. Vale, J. A. Beall, K. D. Irwin, and J. N. 

UUom, Appl. Phys. Lett. 86. 173508 (2005). 

A. M. Savin, M. Prunnila, P. P. Kivinen, J. P. Pekola, J. 

Ahopelto, and A. J. Manninen, Appl. Phys. Lett. 79, 1471 



(2001). 

R. Lconi, G. Arena, M. G. Castcllauo, and G. Torrioli, J. 
Appl. Phys. 85, 3877 (1999). 
"^"^ F. Giazotto, F. Taddei, R. Fazio, and F. Beltram, Appl. 
Phys. Lett. 80, 3784 (2002). 

F. Giazotto, F. Taddei, M. Governale, C. CasteUana, R. 

Fazio, and F. Beltram, Phys. Rev. Lett. 97, 197001 (2006). 

F. Giazotto, F. Taddei, P. D'Amico, R. Fazio, and F. Bel- 
tram, Phys. Rev. B 76, 184518 (2007). 
3° J. P. Pckola, F. Giazotto, and O.-P. Saira, Phys. Rev. Lett. 

98, 037201 (2007). 
31 O.-P.Saira, M. Meschke, F. Giazotto, A. M. Savin, M. 

Mottonen, and J. P. Pekola, Phys. Rev. Lett. 99, 027203 

(2007). 

^'^ J. P. Pckola and F. W. J. Hekking, Phys. Rev. Lett. 98, 

210604 (2007). 

33 A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 

73, 299 (1977) [Sov. Phys. JETP 46, 155 (1977)]. 
3'' N. B.Kopnin, Theory of Nonequilibrium Superconductivity 

(Clarendon, Oxford, 2001). 

35 J. Voutilaincn, T. T. Hcikkila, and N.B. Kopnin, Phys. 
Rev. B 72, 054505 (2005). 

36 R. C. Dynes, J. P. Garno, G. B. Hertel, and T. P. Orlando, 
Phys. Rev. Lett. 53, 2437 (1984). 

3^" J. P. Pekola, T. T. Heikkila, A. M. Savin, J. T. Flyktman, 
F. Giazotto, and F. W. J. Hekking, Phys. Rev. Lett. 92, 
056804 (2004). 



